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Abstract 

Mycobacterium tuberculosis is an important human pathogen, and yet diagnosis remains challenging. Little research has 
focused on the impact of M. tuberculosis on the gut microbiota, despite the significant immunological and homeostatic 
functions of the gastrointestinal tract. To determine the effect of M. tuberculosis infection on the gut microbiota, we 
followed mice from M. tuberculosis aerosol infection until death, using 16S rRNA sequencing. We saw a rapid change in the 
gut microbiota in response to infection, with all mice showing a loss and then recovery of microbial community diversity, 
and found that pre-infection samples clustered separately from post-infection samples, using ecological beta-diversity 
measures. The effect on the fecal microbiota was observed as rapidly as six days following lung infection. Analysis of 
additional mice infected by a different M. tuberculosis strain corroborated these results, together demonstrating that the 
mouse gut microbiota significantly changes with M. tuberculosis infection. 
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Introduction 

Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), 
infects one third of the world's population, and caused 8.7 million 
new cases and 1.4 million deaths in 2010 alone [1]. Although the 
mouse immunopathologic response to TB differs from the 
response in humans, the mouse model is frequendy used to study 
the virulence of different strains of M . tuberculosis, and to assess TB 
drug efficacy [2]. Upon aerosol infection, the bacteria promptly 
replicate in the lungs, and disseminated bacilli are observed in the 
spleen and liver two to four weeks later. Beginning at four weeks 
post-infection, the mouse adaptive immune system achieves partial 
control of the infection, resulting in a plateau in bacterial burden 
[3]. The mouse will develop granuloma-like lesions in the lung 
and, with a high infecting dose, will eventually succumb to the 
bacteria and die. 

One aspect of TB infection that remains largely unexplored is 
the role of the resident microbiota (the microorganisms that 
collectively live on or in mammals). It is now well recognized that 
the gastrointestinal microbiota maintain a complex, reciprocal 
relationship with the host immune system [4,5,6,7]. Moreover, 
differences in fecal microbiota composition and functional 
potential have been identified in individuals with various disease 
states when compared to healthy individuals, including inflamma- 
tory bowel disease, arthritis, type 2 diabetes, and asthma [8,9,10]. 
The gut microbiota also plays a prominent role in both enteric 
bacterial and viral infections [11,12,13,14]. However, there have 



been few studies focused on changes in the gut microbiota in 
response to TB, especially in the absence of antibiotics. 

To date, only one other study has focused on the gut microbiota 
during TB infection. Dubourg and colleagues utilized culture- 
dependent and -independent methods to evaluate bacterial and 
fungal diversity within a single patient at a single time, and found 
that there was an impoverished community dominated by only a 
few phylotypes [15]. However, this patient had been treated with 
broad-spectrum antibiotics for four months, which can gready 
alter the gut microbiota. Another study concluded that there was a 
difference between the sputum microbial composition of TB 
patients and healthy controls [16]. Conversely, a different paper 
focused on human sputum samples found no difference between 
TB patients and healthy controls [17]. 

Given the interaction between the resident microbiota and the 
immune system, as well as competition among the microbiota and 
invading pathogens, we hypothesized that Al. tuberculosis would 
cause a shift in the gut microbiota. Here, we present a longitudinal 
survey using 16S ribosomal gene sequencing of the gut microbiota 
in a mouse model for TB. We assessed bacterial community 
composition and diversity prior to infection with M. tuberculosis 
CDC 1551 and throughout infection until death. Further, we 
evaluated the gut microbiota from additional mice infected by a 
different M. tuberculosis strain (H37Rv) for comparison to the 
longitudinal study. 
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Figure 1. Community structure of individual M. tuberculosis CDC1551 infected mice over time. (A) Survival time in days post-infection for 
each mouse. (B) Phylogenetic profile of bacterial genera. Stacked bar charts in chronological order for each mouse of the 1 8 main genera identified 
based on >1% abundance present in at least two samples. Unclassified sequences are not shown. Black colored bars along x-axis indicate samples 
taken prior to infection, while red colored bars indicate post-infection. Each group represents an individual mouse, followed to death. The mice are 
represented sequentially, with mouse 1 on the left, and mouse 5 on the right. (C) Community diversity in each sample as measured by the Shannon 
diversity index, plotted against the percent survival time. 
doi:1 0.1 371 /journal.pone.0097048.g001 
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Figure 2. Composition of the gut microbiota significantly changes with M. tuberculosis CDC1551 infection. (A) Unweighted and (B) 
weighted Unifrac measures of beta-diversity visualized using Principle Coordinate Analysis (PCoA) following individual mice over time with M. 
tuberculosis CDC1551 infection. Blue dots indicate samples collected pre-infection. Red dots indicate samples collected post-infection. Variance for 
first two component axes is shown as percent of total variance. An analysis of molecular variance (AMOVA) was performed to test whether the 
separation of uninfected and TB-infected samples was statistically significant. In both unweighted and weighted Unifrac measures, there was a 
statistically significant difference (p<0.001). (C) Network analysis of OTUs partitioned among samples, using a five sequence cutoff, and colored by 
phylum. 

doi:1 0.1 371 /journal.pone.0097048.g002 



Results 

Compositional changes in the gut microbiota during M. 
tuberculosis infection 

To study the impact of M. tuberculosis infection on the gut 
microbiota, we infected Balb/c mice with M. tuberculosis 
CDC 1551, and monitored them until death. We collected fecal 
samples at time points prior to infection (pre-infection) and 
throughout infection (post-infection), selecting for analysis three 
pre-infection samples as controls, as well as samples from the first 
two weeks post-infection, the last two weeks prior to death and 
once per month in between (Figure la, Table SI and Figure 
SI). The fecal microbiota was characterized by 454 pyrosequenc- 
ing of bacterial 16S rRNA gene amplicons (VI -V2 region) from 
the five infected mice. A total of 297,156 high-quality sequences 
were generated, corresponding to an average 6,322 reads per 
sample with an average length of 250 base pairs. These sequences 
were clustered into operational taxonomic units (OTUs) at 97% 
pairwise identity and taxonomically classified using the greengenes 
database [18]. The most abundant genera are shown in 
Figure lb. 



We further analyzed overall community diversity using the 
Shannon diversity index, a common ecological diversity measure, 
which takes into account both the number of species (OTUs) 
present and their relative abundance (Figure lc). There was an 
initial decrease in diversity in all mice post-infection, followed by a 
recovery in diversity until death or one week prior to death. This 
was true even for mouse 2, which survived 73 days longer than any 
of the other mice. These trends were also observed using the 
Inverted Simpson diversity index, which takes into account 
community richness, abundance, and is less sensitive to rare 
OTUs compared to the Shannon diversity index (Figure S2). 

Gut community composition and structure differ based 
on infection status 

To identify samples with similar microbial community structure 
and composition, we implemented multidimensional cluster 
analysis based on the weighted and unweighted UniFrac distances 
(Figure 2a and Figure 2b). UniFrac is a phylogenetically-aware 
measure of beta-diversity that can be used to compare OTU 
structure and community diversity. The weighted measure takes 
into account phylogenetic tree branch length. Both measures 
showed clear clustering among the uninfected samples taken pre- 
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Figure 3. Differentially abundant OTUs identified between pre-infection and post-infection. OTUs are ordered by consensus taxonomic 
classification, with OTUs scaled by relative abundances for each row ranging from low relative abundance (blue) to high relative abundance (red). 
doi:1 0.1 371 /journal.pone.0097048.g003 



infection and the infected samples collected post-infection. We 
utilized an analysis of molecular variance (AMOVA), a statistical 
model similar to analysis of variance that is used to analyze 
differences in genetic diversity, to test whether the pre-infection 
and post-infection samples were statistically different, and found 
both weighted and unweighted UniFrac measures were signifi- 
cantly different (/><0.001). We further utilized a network analysis 
to evaluate how the OTUs were partitioned among samples, 
lending additional support for the separation of pre- vs. post- 
infection samples. Strikingly, the first two weeks after infection 
were found to be intermediate between the samples taken before 
infection and later in infection, with the Firmicutes distinguishing 
between the two conditions (Figure 2c). 



We next implemented a categorical analysis to evaluate whether 
specific bacterial OTUs discriminated between pre-infection and 
the infected samples collected post-infection. Eighty-eight OTUs 
were found to be significandy differential (q<0.01); the majority 
belonged to the Firmicutes, specifically within the order Clos- 
tridiales (Figure 3, Table S2) [19]. The bulk of these OTUs 
classified within the Lachnospiraceae family, 3 1 affiliated with the 
unclassified Lachnospiraceae OTU2087 clade, and the Rumino- 
coccaceae family. Interestingly, all of these significantly different 
OTUs were more abundant pre-infection compared to post- 
infection. 
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Figure 4. Phylogenetic profile of bacterial genera for uninfected and M. tuberculosis H37Rv infected mice. Stacked bar charts for 
uninfected and H37Rv-infected mice of the 16 main genera identified based on >1% abundance present in at least two samples. Unclassified 
sequences are not shown. The black colored bar along x-axis indicates the five uninfected mice, while the red colored bar indicates mice infected 
with H37Rv. 

doi:1 0.1 371 /journal.pone.0097048.g004 



Distinct changes in the gut community is independent of 
M. tuberculosis strain 

To determine if our previous results could be replicated with a 
different M. tuberculosis strain, we collected fecal samples from a 
single time point (46 days post infection) from M. tuberculosis 
H37Rv infected Balb/c mice and age-matched uninfected Balb/c 
mice that had been in the same facility for the same amount of 
time. Each group of mice was kept in two different cages to rule 
out variations in caging conditions. As described for the previous 
study, the fecal microbiota was characterized by 454 pyrose- 
quencing of bacterial 16S rRNA gene amplicons (V1-V2 region) 
for this set of mice. A total of 148,466 high-quality sequences were 
generated, corresponding to an average 14,847 reads per sample, 
and were analyzed with the first experiment dataset to maintain 
consistency in OTU clustering (see Methods). The most abundant 
genera are shown in Figure 4. We implemented multidimen- 
sional cluster analysis based on the weighted and unweighted 
UniFrac distances and found distinct clustering between infected 
and uninfected samples, congruent with the results for the 
longitudinal study (Figure 5a, Figure 5b). Network analysis 
similarly confirmed these observations (Figure 5c). These results 
corroborate our earlier observations of distinct clustering between 
infected and uninfected animals. 



Additionally, we found considerable overlap in the phylogenetic 
affiliation of the differentially abundant OTUs, which discrimi- 
nated between uninfected and H37Rv-infected samples, compared 
to the CDC 1551 infection longitudinal study. Seventy-three 
OTUs were found to be differential using the categorical analysis 
(q<0.01), the vast majority similarly within the Lachnospiraceae 
and Ruminococcaceae families (Figure 6, Table S3). Remark- 
ably, the phylogenetic affiliations of the discriminatory OTUs 
mirrored those found in the CDC 1551 infection longitudinal study 
and were highly abundant in the uninfected group, including 
OTU2087 (family Lachnospiraceae), OTU 1995 (order Clostri- 
diales), OTU 1998 (genus Catabacteriaceae), OTU2159 (family 
Ruminococcaceae), OTU2166 (genus Clostridium), and OTU2176 
(genus Oscillospira). However, while there was high congruence in 
the phylogenetic lineages observed between the two M. tuberculosis 
studies, only five OTUs specifically overlapped in both exper- 
iments (Table S2 and Table S3). 

Discussion 

We have monitored changes in the composition of the mouse 
gut microbiota from pre-infection to death, and shown that there is 
a clear difference between the microbial communities of infected 
and uninfected mice, results that have been confirmed by the use 
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Figure 5. Gut microbiota composition of M. tuberculosis H37Rv infected mice is significantly different from uninfected mice. (A) 

Unweighted and (B) weighted Unifrac measures of beta-diversity visualized using Principle Coordinate Analysis (PCoA) for the comparison of H37Rv- 
infected mice to uninfected mice at a single time point. Blue dots indicate samples collected pre-infection. Red dots indicate samples collected post- 
infection. Variance for first two component axes is shown as percent of total variance. In both unweighted and weighted Unifrac measures, there was 
a statistically significant difference (AMOVA p<0.005). (C) Network analysis of OTUs partitioned among samples, using a five sequence cutoff, and 
colored by phylum. 
doi:1 0.1 371 /journal.pone.0097048.g005 



of two different strains of M. tuberculosis. These changes occur by 
day six post aerosol infection, indicating that this shift is very rapid. 
The prompt shift in the community after aerosol infection suggests 
that the gut microbiota is modulated during TB infection and may 
be responding to host immunological changes [20]. Furthermore, 
there is a consistent trend in these alterations, with a decrease in 
overall community diversity and discrete shifts in specific microbial 
taxa, despite the difference in survival time for the five mice 
studied. This trend occurred even in mouse 2, which survived 73 
days longer than the other mice. Overall, we identified a number 
of OTUs significantly different between infected and uninfected 
mice, regardless of the infecting strain of M. tuberculosis. Mycobac- 
terial DNA was not detected in any of these samples with the 16S 
rRNA pyrosequencing, but PGR of the IS6110 insertion sequence, 
which is specific for M. tuberculosis, in the last sample collected prior 
to death identified TB in only three of the five mice (data not 
shown). We posit that the observed microbiota changes are not be 
mediated directly by the presence of M. tuberculosis in the gut, but 
instead represent crosstalk between the resident microbiota and 
the mucosal immune system. In fact, the minimum in diversity 
approximately corresponds to the time when the adaptive immune 
system begins to achieve a plateau in bacterial burden [21]. Thus, 



the loss of diversity could be a result of immune system activation, 
with a recovery once bacterial burden and immune activity have 
reached equilibrium. 

We observed significant differences in the relative abundance of 
members within the Lachnospiraceae and Ruminococcaceae 
families (Clostridiales) in the two studies, with higher relative 
abundances in both the uninfected and pre-infected mice. It is 
known that some members of clusters IV and XlVa of the genus 
Clostridium (a genus in the Clostridiales) induce regulatory T cells 
[22]. Members of the Lachnospiracaceae have also been shown to 
decrease in abundance in inflammatory bowel disease (IBD) and 
colitis [23,24] . Furthermore, a recent study found that a decrease 
in abundance of OTU 2087 was correlated with exacerbated 
asthma [25]. Future work to delineate immune-modulatory 
members of the Clostridiales, particularly those members associ- 
ated with regulatory T cells and which have a role in M. tuberculosis 
immunity, will be critical. 

While not as prominent as the differential Clostridiales OTUs, we 
additionally observed a subset of Bacteroidales OTUs differentially 
abundant in pre-infected and uninfected samples in the two 
studies. Bacteroides spp. have been shown to play important anti- 
inflammatory roles in inhibiting activation of the NF-kB pathway 
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Figure 6. Differentially abundant OTUs identified between uninfected and M. tuberculosis H37Rv infected mice. OTUs are ordered by 
consensus taxonomic classification, with OTUs scaled by relative abundances for each row ranging from low relative abundance (blue) to high 
relative abundance (red). 
doi:1 0.1 371 /journal.pone.0097048.g006 



and induction of IL- 1 0-producing T cells [6]. Furthermore, 
Bacteroides fragilis has been shown to modulate the T-helper type 1 / 
2 (Thl/Th2) balance, another aspect of the immune system 
critical for the control of M. tuberculosis [26]. Further studies to 
characterize the relationship of specific Bacteroides species to the 
immune system are clearly needed. 

Together, our results establish that the gut microbiota of mice 
changes significantly following aerosol infection by M. tuberculosis, 
and these differences may be related to the immune signaling from 
lung to gut. These changes begin as early as six days post infection 
and are characterized by an initial loss in diversity, which recovers 
to a significantly different composition. During this time, there are 
many alterations in the relative abundance of a number of OTUs, 
most significantly a decrease in members of the Clostridiales and 
Bacteroidales. These observations have important implications for 
our understanding of the interplay between the immune system 
and the gut microbiota. Further investigation into the interaction 
of the host immune system and the gut microbiota could lead to a 
more mechanistic understanding of gut immune functioning and 
the role the microbial fraction plays during TB infection. 

Materials and Methods 

Ethics Statement 

This study was carried out in strict accordance with the 
recommendations in the Guide for the Care and Use of 



Laboratory Animals of the National Institutes of Health. The 
Animal Protocol MO11M120 was approved by the Institutional 
Animal Care and Use Committee of the Johns Hopkins 
University, and covers all animal procedures, including time-to- 
death. All efforts were made to minimize suffering. 

Bacterial Strains 

Mycobacterium tuberculosis H37Rv and CDC1551 were routinely 
grown in Middlebrook 7H9 broth (Fisher Scientific, Waltham, 
MA) supplemented with 0.5% glycerol (Sigma, St. Louis, MO), 
0.05% Tween 80 (Sigma) and 10% oleic acid-albumin-dextrose- 
catalase (OADC; Fisher Scientific) at 37°C with agitation. 

Animals 

For the M. tuberculosis CDC 1551 experiment, 4-6 week old 
female Balb/c mice were purchased from Charles River 
(Wilmington, MA). The five mice used for stool collection were 
toe tattooed for identification purposes and housed in the same 
cage. Mice were monitored daily as part of a time-to-death 
experiment and were provided with gel cups broken up feed 
biscuits when they became cachectic. For the M. tuberculosis H37Rv 
experiment, 19-20 gram female Balb/ c mice were purchased from 
Charles River. Mice were housed in cages containing 5 animals. 
For each group, 3 mice were randomly selected from one cage, 
and 2 from another, for stool collection. For colony forming unit 
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(CFU) counts of lungs and spleens, mice were euthanized by 
cervical dislocation. 

Stool Collection and Storage 

For stool collection, each mouse was temporarily placed alone 
in a clean container until there was approximately 0.15 g stool in 
the container. The mouse was removed and the stool collected and 
stored in O-ring sealed tubes (Simport, Beloeil, Canada) at 
-80°C. 

Aerosol Infection 

Mice were infected with log-phase broth cultures diluted in 
phosphate-buffered saline (PBS) using the Middlebrook inhalation 
exposure system (Glas-Col, Terre Haute, IN) in a single run. 
Bacterial burden was determined by sacrificing 3 mice per time 
point and enumerating CFUs on selective 7H11 plates (VWR, 
Radnor, PA). 

DNA Extraction and 16S rRNA Sequencing 

Total DNA was extracted from 0.15 grams of stool using the ZR 
Fecal DNA isolation kit (ZYMO Research Corp., Irvine, CA) with 
modifications including an enzymatic pre-treatment with muta- 
nolysin (5000 units ml -1 , Sigma) and lysozyme (100 mg ml -1 , 
Sigma) in conjunction with aggressive bead beating using 0.1 mm 
glass beads (BioSpec, Bartlesville, OK) and a bead-beater 
(BioSpec). Barcoded primers [27] were used to amplify the 
bacterial 16S rRNA gene region 27F-338R from 50 ng of purified 
DNA using AccuPrime High Fidelity DNA polymerase (Invitro- 
gen, Grand Island, NY) in a total reaction volume of 25 uL. 
Reactions were run in a PTC- 100 thermal controller (MJ 
Research, Waltham, MA) using the following cycling parameters: 
5 min of denaturation at 95°C, followed by 25 cycles of 30 s at 
95°C, 30 s at 55°C, and 60 s at 68°C, with a final extension at 
72°C for 7 min. Amplicons were quantified using the Quant-iT 
PicoGreen dsDNA assay and equimolar amounts (100 ng) of the 
PGR product were mixed in a single tube. Amplification primers 
and reaction buffer were removed from each sample using the 
AMPure Kit (Agencourt, Brea, CA). The purified amplicon 
mixtures were sequenced by 454 FLX pyrosequencing using 454 
Life Sciences primer A by the Genomics Resource Center at the 
Institute for Genome Sciences, University of Maryland School of 
Medicine. 

16S rRNA sequence processing and analysis 

Sequences generated from pyrosequencing of bacterial 16S 
rRNA gene amplicons from all mouse experimental groups were 
processed using mothur, version 1.27, according to the standard 
pipeline oudined in [28,29]. Sequences were denoised, trimmed, 
quality and chimera checked using de novo uchime, and clustered 
into operational taxonomic units (OTUs) at 97% pairwise identity 
and aligned to the Silva reference alignment [30]. Taxonomic 
classifications were assigned using the naive Bayesian classifier 
with the May 2011 release of the greengenes database [18]. 
Rarefied OTUs (randomly subsampled to normalize sequence 
counts) were used to calculate community diversity for each 
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